THE SIMPLEST MODEL OF SPATIALLY DISTRIBUTED POPULATION WITH 
REASONABLE MIGRATION OF ORGANISMS 
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The simplest model of a smart spatial redistribution of individuals is proposed. A single-species 
population is considered, to be composed of two discrete subpopulations inhabiting two stations; 
migration is a transfer between them. The migration is not random and yields the maximization of 
net reproduction, with respect to the transaction costs. The organisms are supposed to be globally 
informed. Discrete time model is studied, since it shows all the features of a smart migrations, while 
the continuous time case brings no new knowledge but the technical problems. Some properties of 
the model are studied and discussed. 

PACS numbers: 87.23.Cc, 87.23.Kg, 87.23.-n 



I. INTRODUCTION 

Modelling of the dynamics of biological takes start 
from the works by V.Volterra jME^I (see also The 
approach follows, in brief, a chemical kinetics, where re- 
production, death and other types of interaction of organ- 
isms are described with law of mass action. A spatially 
distributed population is modelled with the "reaction - 
diffusion" system, where diffusion is supposed to describe 
the redistribution of organisms over a space. Good co- 
incidence of the solutions of such "reaction - diffusion" 
systems to observed dynamics of populations and com- 
munities just masks the serious problem. 

The point is that the diversity and abundance of possi- 
ble (structurally stable) regimes of such dynamic systems 
exceed drastically any really available family of trajecto- 
ries of real systems. Thus, one always is able to match 
a differential equation (or a system of differential equa- 
tion, either ordinary, or partial differential equations) 
to the dynamic behavior of any really observed system. 
Doubtlessly, living organisms, even microorganisms, dif- 
fer drastically in their "microscopic" behavior from the 
chemical substances and relevant chemical reactions. 

Vito Volterra, the founder of the mathematical biol- 
ogy of populations understood pretty well this discrep- 
ancy. Later, the comprehension of limitations of this 
chemical methodology fell off. The situation is going 
worse when one tends to model a dynamics of spa- 
tially distributed populations and communities. Ba- 
sically, "reaction - diffusion" systems make the ba- 
sis for modelling of such spatiall y d istributed entities 

[1 II II II II ES 111 111 111 El 11 113 Famous so- 
viet mathematician Andrew Kolmogorov studied in de- 
tail such type of equations ^3 ■ 

A methodology of modelling of spatially distributed 
populations based of the reaction - diffusion system has 
the great discrepancy. It constrains significantly the indi- 
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vidual (behavioral) properties of organisms under consid- 
eration: one must suppose that the organisms move over 
a space randomly and aimlessly . ObviouslVjSuch as- 
sumption does not hold true (see, e.g. HI^E3)j even 
for microorganisms 0, HI H^, 0| • The assumption to- 
wards the idle transfers of organisms in space is obviously 
less favorable for the higher organised species. 

Modelling methodology based on evolution principles 
is the way to pass over the discrepancy mentioned above. 
This approach takes the origin from the evolution studies 
of J.B.S.Haldane HI]- This is the most general principle 
prescribing the way biological systems evolve. In brief, 
it force to evolve a system toward the maximization of 
net reproduction. This latter is an average number of 
per capita descendants survived at the course of a series 
of reproductions over an arbitrary long generation line 
[T^ IT3 . 113- Later, they found this principle to be 
even more general, than just a biological one. Indeed, the 
principle holds true for a ny sy stem where the inheritance 
takes place H E S HI 111113 ■ 

A consistent and rational implementation of this prin- 
ciple faces the problem of a lack of knowledge of how the 
specific biological issues impact the survival of a species. 
In turn, the question arises, what is an entity to be evolv- 
ing? An ordinary answer on this question is that the 
species is an evolving entity. Actually, the situation looks 
more complicated; not discussing this problem in detail, 
further we shall follow this idea. The principle formulated 
above yields the following rule for the model implemented 
below: evolution optimality in spatial distribution of or- 
ganisms is equivalent to the maximization of (an average) 
net reproduction over space, with respect to the evolution 
trade off for such redistribution. Some further details on 
this issue could be found in HI QUI H^. 



II. BASIC MODEL OF THE SMART 
MIGRATION 

Consider a population inhabiting two stations; hence, 
the population consists of two subpopulations. Any 
movements of individuals within a station must be ne- 
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glected; thus, only the transfer from station to station is 
considered as a migration act. No spatial effects in the 
population dynamics are presumed, for each subpopula- 
tion, as soon as no migration occurs. We shall consider a 
discrete time model; continuous modelling is possible, as 
well, but it brings no new issues but the serious technical 
difhculties. 

Further, the dynamics of each subpopulation is sup- 
posed to follow the Verchulst's equation IsllslliollilLlisll. 

Namely, let Nt and denote the abundance of the first 
subpopulation (of the second one, respectively) , so that 



and 



Nt+i ^a-Nt-b-N^ 



Mt+i = c • Mt ~ d ■ Af/ 



(la) 



(lb) 



respectively, for the case of the absence of migration. We 
shall implement the migration into consideration later. 
Here a and c represent fertility of the relevant subpop- 
ulation, while b and d show the effect of density depen- 
dent competition within a subpopulation, each. It must 
be said, that unlike for the Verchulst's equation of a sin- 
gle population, a student may not change the equations 
(|lal llb|l for the dimensionless form, eliminating the co- 
efficients b and d. The point is that the migration effects 
would break down such transformation. The functions 



kr{Nt\ 



b-Nt 



and 



ki{Mt) = c ~ d ■ Mt 



(2a) 



(2b) 



are the reproduction coefficients, in relevant stations, re- 
spectively. 

Migration is a transfer of a part of subpopulation from 
one station into the other. Migration itself affects the 
reproduction, survival and other vital functions of an or- 
ganism. All these issues will be integrated into the pa- 
rameter called the cost of migration /i: no negative im- 
pact on the reproduction, survival and other vital func- 
tions of an organism is observed, as /i = 0. Otherwise, 
the growth of fi yields a decay in survival of organisms. 
Rather simple way to introduce the migration cost is pro- 
posed in [m liglls^ : 

P = exp{- ^} , 

so that p, < p < I becomes a probability of the trans- 
fer from station to station with no damages for further 
reproduction. 

Migration starts, as soon as the conditions of life (mea- 
sured in the units of reproduction coefficient ^) "there" 
becomes better, than "here" . This narrative issue could 
be easily transformed into the formal way. Migration 
from TV-station (from M-station, respectively) starts, if 
p-ki{M) > kr{N) (if ki{M) <p-kr{N), respectively) for 



the given time moment t, t = 1,2,3,.... We shall sup- 
pose, the migration runs as a single act, entirely, and the 
reproduction always follows the migration. Thus, a life 
cycle of a population consists of two stages: the former 
is a migration (if it takes place at the given time moment 
t), and the latter is the reproduction ran according to 
([T|). with the abundances appeared due to migration, at 
each station, independently. 

Schematically, the model works in three steps. 
The first step consists in the determination of the fact 
of migration, i.e. in the comparison of reproduction co- 
efficients "there" and "here": 

• if p ■ ki{M) > kr{N), then the migration runs from 
A^-station to Af-station; 

• if ki{M) < p ■ kr{N), then the migration runs from 
Af-station to A^-station; 

• otherwise no migration takes place. 

The second step consists in the determination of the 
migration flux. Migration yields the change of life con- 
ditions measured in the units of reproduction coefficient 
(0). An emigration results in the growth of the coeffi- 
cient, reciprocally, immigration results in the decrease of 
that latter, since the coefficients ^ are supposed to be a 
linear descending function of the abundance A^ (or Af). 
The number of migrating individuals tends to equalize 
the reproductions coefficients, due to an abundance re- 
distribution: 

ki{M - A) =p-kriN+p-A) 

for the case of migration from Af -station to A^-station, 
and vice versa: 

p-ki{M + p-A) =kriN -A) . 

Here the term p ■ A represents the fact of mortality of 
individuals at the course of migration; this is the way to 
account the transfer cost, for this model. The migration 
flux, then, is determined according to 



pc-a + bN ~ pdM 



b + p^d 

for the migration from A^-station to Af-station, and 
pa - c + dM - pbN 



Anm — 



d + p% 



(3a) 



(3b) 



for the inverse migration. 

Finally, the third step consists in the reproduction of 
organisms in both subpopulations independently, with 
respect to the abundance (iVj, or A^/t, respectively) re- 
sulted from the migration: 



t+i 



and 



t+i 



a-Nf-c-Nf 



c-Mt-d- Mf 



(4a) 



(4b) 
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It should be stressed, that a reproduction in a station 
(i.e., at the subpopulations) runs independently, with re- 
spect to the abundances occurred due to the migration. 
If no migration takes place for some time moment t, then 
Nt = Nt (Mt = Mt, respectively). 



III. SOME PROPERTIES OF THE MODEL 

The model (^-0J of the smart migration exhibits vari- 
ous dynamic properties, not observed at the similar pop- 
ulation dynamics models with no migration. The dynam- 
ics of two (independent) subpopulations runs inside the 
rectangular [0,a/6] x [0, c/c?], if no migration occurs. As 
soon as p > 0, the dynamics runs at the dovetail shown 
in FigH If a trajectory (i.e. a point representing a cou- 
ple of abundances {Nt,Mt)) reaches the area 1, then no 
migration occurs at that time moment and no migra- 
tion will take place, while the trajectory remains at this 
area. The area is cut off with two solid bold lines. These 
are determined by the equations 

Aa/ jv = and Is.mm = . 

Areas 2 and 3 are the space dynamically available by each 
subpopulation, independently, when no migration takes 
place. 

On the contrary, the migration expands the dynam- 
ically reachable area. Such expansion results from the 
smart migration: individuals emigrate from the overpop- 
ulated station improving the survival of the entire popu- 
lation. Areas 2 and 5 in Fig^show this expansion. For 
any point (iVj, Mt) G Jl, where is a union of these two 
areas, the migration gives the projection of this point 
on the bold solid line bordering the area 1, parallel to 
the solid dashed line shown in this Figure. Thin arrow 
in this Figure illustrates such projection. Similarly, if a 
point [Nt, Mt) S $ belongs to the union $ of areas 3 and 
4, then the migration maps it on the (lower) solid bold 
line bordering area 1. Similarly, it is executed in parallel 
to the thin arrow, shown in the figure. 

The tangent of the dashed line bordering the area 5 
is equal —p~^\ similarly, for the area 4 the tangent is 
equal to —p. The area 1 expands, as p — > occupying 
the entire rectangular [0, a/5] x [0, c/d], for p = Q. The 
areas 4 and 5 become a (semi-infinite) strip each, of the 
permanent width. This fact differs the situation of the 
complete absence of a migration from the infinitely big 
migration cost /x. This area collapses into the line (into 
the intercept, to be exact) defined the the equation 

hN-dM + a-c = {), 

as p 1. 

Next, the migration expands the allowed parameter 
values. The parameters a and c must meet the constraint 

l<a<4 and l<c<4. 



in migration free models H IH IH IH HI El [i^ . The 
equation Ijlal) (or (jlbp . respectively) exhibits elimination 
of a population, as a (c, respectively) exceeds 4: maxi- 
mum of aNt—bNf may be greater than a/b. The non-zero 
migration expands the range of the parameter, provided 
that the overpopulation in one station will be compen- 
sated by the emigration into another. 

Equation (|la|l (or (|lbp , respectively) shows various dy- 
namic patterns, for various parameter a (or parameter c, 
respectively) value. The diversity of limit regimes of 
varies from a stable fixed point to a strange attractor as 
a limit manifold. Similar diversity of regimes could be 
found at the equation system Ijlal llb|l , when the migra- 
tion occurs. In addition, the system (^-0)) exhibits some 
other regimes, that are not possible with no migration. 

Migration provides a redistribution of individuals over 
a space. The redistribution is not random; it results in 
maximization of the average (over two stations) repro- 
duction coefhcient kr{N) + ki{M). The reproduction co- 
efficient reaches the maximum in one step; this point 
results both from the global information accessibility of 
the data concerning the environmental conditions (coeffi- 
cients a, b, c and d), population density (these are Nt and 
Mt in both stations, at every time moment t), and the 
transfer cost fj, (or the probability of the successful trans- 
fer p). If < p < 1, the migration results a decrease of a 
general abundance, since a part of individuals (namely, 
p • A) is eliminated, at each time step t. Migration yields 
no elimination of individuals, as p = 1. 

Various limit regimes occurred due to smart migration 
effect may be observed. There exists the stable perma- 
nent one-side directed flux of individuals. For example, 
if c is big enough, and simultaneously b is sufficiently 
small (thus increasing the environmental capacity of A'^- 
station) , one may observe the permanent one-way migra- 
tion flux limit regime. Indeed, such regime is determined 
by the equation 

f A^* = a- (A^* +pA*) - 6- (A^* -l-pA*)^ , , 
\ M* = c • {M* - A*) - d • (M* - A*f , ^ ' 

with A* determined by l|3b)l : 

_ pa - c + dM* - pbN* 
^ d + p^b 

To figure out the impact of the smart migration on the 
population dynamics, let's consider a particular case of 
b — d — 1^ and p = 1. Such choice of the parameters 
means that the subpopulations differ in their growth rate, 
only, and no losses of an abundance are resulted from a 
migration. The equality p = 1 also means a supreme 
mobility of an individual. Here A becomes equal to 

A = a- c-HM-iV, 

and the system (0 changes for 

[ N ^ a-{M + X)-{M + \f , . 

\ M ^ c- (N - X) - {N - Xf , ^ ' 
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FIG. 1: Space available for the dynamic trajectories of the model, in case of non-zero migration; 1 - no migration area, 2, 5 - 
areas of migration from A^'-station to M-station, 3,4 ~ areas of migration from -station to A^'-station. Horizontal axis shows 
the abundance at A^-station, vertical one shows similar abundance in M-station. 



with A = a — c. Here the asterisks are omitted, since 
it makes no confusion. This is the system of two poly- 
nomials of two variables (Af, iV) of (formally) power 2 



J 



each. To solve a system of two polynomials of two vari- 
ables, one must develop the resultant of the system. This 
former is the determinant of the following matrix: 



/O 1 



1 -(2A + c) 

Vo 1 



[a(Af + A) - {M + A)2] 
1 

A2 + cA + M 
-(2A + c) 



\ 

\a{M + A) - {M + A)2] 



A2 + cA + M ) 



(7) 



r 



for the case Some roots of ((Jj) are the solution of lO. 
It should be stressed, that the first term at the first line 
of the matrix Q becomes non-zero, for the case h ^ d. 

Suppose, (iV*, M*) is the solution of ®, and N* > 0, 
M* > 0. Here the question arises, whether this solution 
is stable. The answer on this question could be obtained 
due to linear approximation analysis. Let Nt — N + i^t 
and Mt — M + fit, where are the small correc- 

tions; substituting such {Nt, Mt) into lO and omitting 
the terms of the second and higher orders, one gets the 
following matrix for linear approximation: 



-2 • (A/ + A) 

(iV - A) 



(8) 



The eigenvalues of matrix (jSJ are less, than 1, when 
\{M + X) ■ {N - X)\ < 1/4. 

Similarly, numerous other limit regimes of different 
complexity and structure could be found, as well. One 



hardly can figure them out explicitly; moreover, there is 
no much sense in detailed determination of these com- 
plex regimes. Computer simulation makes them rather 
obvious. 

Here we present the simplest model of a smart mi- 
gration, strongly opposing the "reaction - diffusion" 
methodology. The model is rather clear and apparent. 
The simplicity results from the peculiar feature of the 
model; this is the case of globally informed individuals. 
Indeed, an individual makes a decision whether it would 
migrate, or would stay in the station, referring to the 
available information concerning the life conditions. So, 
the key issue is what information towards that matter is 
available. 

Let's concentrate on some mathematical issues fol- 
lowed from the smart migration. The first one is that 
the dynamics of a population with smart migration is ir- 
reversible in time. Indeed, the migration at the model 
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(^-"0J is a projection, from mathematical point of view; 
thus, a set of different states are transformed into unique 
one, and there is no regular way to figure out which one 
was preceding the observed abundance resulted from the 
migration act. 

Another important issue is that it expands both per- 
missible phase space, and the parameters values, in com- 
parison to similar models with no migration. It should be 
said, that such effects, probably, also could be observed 
(at least, for some peculiar combinations of the parame- 
ters) for the systems with random, or aimless migration. 
We doubt that the direct comparison of the areas of per- 
missible phase variables, or the parameters, for the case 
of smart vs. the aimless migration makes much sense. 
Obviously, the model is rather simple and specific to pre- 
tend to describe properly any real biological system. The 
specificity here manifests in the global information ac- 
cessibility; in occupation not more than two stations; in 
the absence of "microscopic" consideration of dynamics 
within a station. 

The main purpose of this paper is to show the method- 
ology of the modelling of spatially distributed popula- 
tions (and other biological communities) with no chemi- 
cal analogies, rather than to find out some peculiar dy- 
namic regime pretending to match exactly a dynamics of 
a real biological community. 

The model provides that the life conditions at the resi- 
dence station are converted into a couple of parameters a 
and h (c and d, respectively). Besides, it is supposed that 
the density of the subpopulation (or its total abundance 
iV, or M, respectively) is known, as well. Not discussing 
at the moment the details of the detection of population 
density, or other conditional parameters per se, suppose 
that the environmental abundance (density, indeed) and 
other parameters are detectable for an individual. Be- 
sides, the model (P - 0J suggests that similar parame- 
ters, and the abundance are known, at the distant sta- 
tion. Such suggestion makes the individuals to be glob- 
ally informed. A feasibility of such presumption is doubt- 
ful, nevertheless, this assumption is a comm on p lace for 
mathematical population biology [H H 11121 Im I2ll3^ . 

The methodology presented above is rather powerful, 
and provides a researcher with the tool for studying spa- 
tially distributed populations with no artificial and ab- 
solutely unrealistic hypotheses towards the microscopic 
behavior of individuals, i.e. the randomness and aimless- 
ness of their transfers over a space. One sees the follow- 
ing furthering of the approach described here. First, a 
two-species (or several species) communities could be de- 
scribed within the framework of the methodology. Again, 
one should consider a two-station model, where each 
species (say, predators and preys) migrate from station to 
station and back. A dynamics within a station might be 
modelled with the the most common equation (say, with 
Lotka-Volterra equation), thus explicating the effect of 
smart migration in the dynamics of a multi-species com- 
munity. 

Both single-species, and multi-species models with 



smart migration may incorporate various patterns of in- 
formation accessibility, for individuals. The model pre- 
sented above is based on the hypothesis of the global 
accessibility of information to an individual. It means, 
that an individual knows the conditions of life (expressed 
in the coefficients, at the case of the model (n]-0J, both 
at the station of residence, and the station of immigra- 
tion; it knows the density of each subpopulation. Finally, 
the an individual knows the transfer cost, in this case. 

The hypothesis of total lack of information available 
to individuals opposes the idea mentioned above. Here 
an individual operates with the inner, extremely local in- 
formation, when makes a decision on the change of the 
location for another one. In brief, such situation could be 
described like a threshold migration, where the transfer 
act takes place only when the local conditions become 
worse than some individually defined level. Still, the sit- 
uation of the total lack of information is not equivalent to 
random and aimless migration. The difference becomes 
clear, if one considers the situation of the transfer act oc- 
currence, while the local density at the occupancy place 
is still very low. The smart migration under the total 
lock of information would start up, while chemically-like 
diffusive migration will not take place, in such case. 

All these assumptions seem to be too strong and spe- 
cific. In general, the individuals operate with a part 
of information. There are several problems here, both 
of mathematical origin, and of biological essence. The 
first one consists in exact and comprehensive definition 
of what exactly is known to individuals. Next one is the 
discretion between the behavioral patterns supported by 
the reasonable choice of the way to behave, and those 
determined instinctively. Consider the seasonal bird mi- 
gration. Surely, the fact that some species change a re- 
production site for a winter spending site, falls beyond 
the will of a bird, it does not make a matter of reason- 
able of self-made choice: that is the instinct forcing them 
to fiy away. On the other hand, the choice of a pecu- 
liar site to spend a night (if any) is made by the birds 
in a flock reasonably, with respect to the detail features 
and circumstances of the current situation. Finally, the 
problem arises when one tends to determine how far (in 
space) the individuals are able to collect and process the 
information (concerning the living conditions "there"). 

The model - 0J implies that the transfer cost /i 
is symmetrical, and does not depend on the direction of 
migration. It might be so, while more realistic idea is that 
the transfer cost should be unsymmetrical. Evidently, 
the simplest way to figure out the transfer cost is to split 
it into three parts: 

M = Mout (A)+/iin(B)+^(A,B) . (9) 

Here /iout(A) is the transfer cost of successful emigra- 
tion from the station A; /iin(B) is the transfer cost of 
successful introduction into the station B, respectively; 
finally, (5(A, B) is the pure transfer cost from station A 
to station B. Obviously, one should expect to face the 
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asymmetry 

Aiout(A) ^ A*out(B) and Afin(A) ^ A*in(B) , 

in general. Also, a symmetry of pure transfer cost 
(5(A, B) = 5(B, A) is doubtful. It should be said, that 
this point makes no problem in its implementation at 
the model JTl-|i|). 

Another significant constraint of the model (^-0J is 
the spatial structure limited with two stations. Indeed, 
an expansion of the approach presented above for the 
case of several stations, and, ultimately, for a contin- 
uous, or quasi-continuous case of a habitat is strongly 
desirable. Suppose, a population inhabits three stations; 
here we presume the global information accessibility, as 
well. Suppose, further, the conditions (i.e., the abun- 
dances and the parameters) make the situation when the 
individuals from the station A must migrate either to 
station B, or to station C. No one knows exactly, in ad- 
vance, what is a proportion of individuals immigrating 
into the station B vs. those immigrating the station C. 
This is the main obstacle here. There exists the approach 
withdrawing this discrepancy; it is based on the interval 
mathematics 21, 23, 24, 34, 35j. The detailed discussion 
of that issue falls beyond the scope of this paper. 



IV. CONCLUSION 

The model described above implements the method- 
ology of evolution optimality into the problem of the 
modelling of spatially distributed populations. Migration 
causes the growth on net reproduction (which is a repro- 
duction rate, in our case), in average, over the space. The 
model comprises the simplest case of two stations, where 
the spatial distribution is restricted to a transfer of indi- 
viduals from station to station and back; the transfer cost 
is supposed to be symmetrical one. The model shows the 
expansion of the environment capacity, in comparison to 
the case of the migration absence. 
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